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Abstract 

Let T ab = T ba = be a system of differential equations for the components of a 
metric tensor on M m . Suppose that T ab transforms tensorially under the action 
of the diffeomorphism group on metrics and that the covariant divergence of T ab 
vanishes. We then prove that T ab = E ab (L) is the Euler-La grange expression some 
Lagrangian density L provided that T ab is of third order. Our result extends the 
classical works of Cartan, Weyl, Vermeil, Lovelock, and Takens on identifying 
field equations for the metric tensor with the symmetries and conservation laws 
of the Einstein equations. 

Keywords: metric field theories, symmetries of differential equations, conserva- 
tion laws, inverse problem of the calculus of variations. 
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1 Introduction and main results 

The time-honored Noether's theorems [23] establish a correspondence between symme- 
tries and conservation laws for the Euler-Lagrange equations of a classical variational 
problem. Noether's first theorem states that every infinitesimal symmetry of the varia- 
tional problem determines a differential conservation law for the Euler-Lagrange equa- 
tions, and conversely that, under certain mild non- degeneracy conditions, conservation 
laws can be associated with symmetries of the variational problem. Noether's second 
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theorem, in turn, asserts that infinite dimensional symmetry pseudogroups of the varia- 
tional problem involving arbitrary functions of all the independent variables correspond 
to differential identities among the equations. 

In 1977, F. Takens [27J considered the following novel and very distinct aspect of the 
interplay between symmetries, conservation laws and variational principles for systems 
of differential equations: 

Let r be a Lie algebra of vector fields defined on the space of independent 
and dependent variables, and suppose that a system of differential equations 
is invariant under T and that each element in T generates a conservation law 
for the system. Does it then follow that the system arises from a variational 
principle? 

In his original paper Takens studied the question for second order scalar equations, 
systems of linear equations, and equations arising in metric field theories. Takens' 
results for second order scalar equations and for systems of linear equations were sub- 
sequently generalized by Anderson and Pohjanpelto [7], [8], [21]. We refer to [7] in 
particular for more background material and motivation on Takens' problem. 

In addition to the papers listed above, there is an extensive literature dealing with 
the existence of variational principles for systems of differential equations admitting a 
Lie algebra of symmetries and the corresponding conservation laws within the context of 
Noether's second theorem, that is, where the symmetry group is the infinite dimensional 
group of coordinate transformations of the underlying manifold and the conservation 
laws express the vanishing of the covariant divergence (or some variant of it) of the field 
equations. This work is largely directed towards the axiomatic characterization of the 
Einstein equations. The original classification results of Cartan [11] , Vermeil [28], and 
Weyl [29] establish that second order quasi-linear field equations for the metric tensor 
possessing the symmetries and conservation laws of the Einstein equations necessarily 
arise from a variational principle. These results were later extended to general, fully 
nonlinear second order systems for the metric tensor and to third order systems in the 
3-dimensional case by Lovelock, JT6J [TTJ [19], again by a detailed classification of all 
equations with the required properties. A direct proof for second order systems based 
on the analysis of the Helmholtz conditions [2J, i.e., the integrability conditions for 
the existence of a variational principle, can be found in [27]. Lovelock's results were 
later generalized to metric-scalar [13], [13], metric- vector |18] . and metric-bivector [20J 
theories. Anderson [3] subsequently provides a general characterization of second order 
divergence free systems for the metric tensor and an auxiliary independent tensor field, 
subsuming in part the above-mentioned works on combined metric field theories. 

In [13], Horndeski attempts to extend Lovelock's work to general third order equa- 
tions with the symmetries and conservation laws of the Einstein equations. However, 
the treatment for the existence of a variational principle falls short of a comprehen- 
sive result due to a restrictive extraneous assumption that the zeroth order Helmoltz 
conditions for the system be invariant under the action of the diffeomorphism group. 
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The next important step extending Lovelock's work was the introduction of the 
generalized Cotton tensors in [Tj . These tensors and their construction is placed into the 
proper differential geometric context by Anderson p)j as the Euler-Lagrange expressions 
of Lagrangians derived from the Chern-Simons forms of Riemannian geometry. As is 
discovered in [5] , these analogs of the classical Cotton tensors of 3-dimensional conformal 
geometry play a key role in the equivariant inverse problem of the calculus of variations 
for metrics. 

The aim of this paper is to extend the results of Cartan, Vermeil, Weyl, and Lovelock 
to general third order field equations for the metric tensor with the symmetries and 
conservation laws of the Einstein equations. As is intimated by the intricacy of the 
constructions in [TJ, [5], a direct classification of third order equations along the lines 
of the original works on the subject would be a formidable undertaking. However, as 
is well known, the existence of a local Lagrangian for a system of differential equations 
is ensured by the vanishing of the classical Helmholtz conditions for the equations, and 
our problem is rendered tractable by an analysis of these conditions for systems sharing 
the properties of the Einstein equations. 

To describe our results more precisely, write (x a ) for the coordinates on M. m . A 
metric 

g = g a bdx a <g> dx b 

on M. m is a symmetric type (0,2) tensor field with g = det(g a b) + 0, where the g a b 
stand for the components of g. In this paper we consider metrics of fixed but arbitrary 
signature. The action of the Lie algebra X(M. m ) of vector fields on metric tensors via 
pull-back gives rise to the infinitesimal transformation group 

o-^-^-^i;^^^} (L1) 

on the coordinate space G = {(x l , g& c )}. 

The metric tensor g is subject to a system of kth order partial differential equations 

T — T^ ^ (x , g c d, gcd,i\ j gcd,i\i 2 i ■ • ■ t gcd,i\i2- -ik ) — ^' ft, 6 — 1, . . . , 771, 

where g c d,hi 2 -ii denotes the derivative of g c d with respect to the independent variables 
l . The operator T ab is locally variational if it can be written in some 
neighborhood of each point of its domain as the Euler-Lagrange expression 

T - = E"\L) - - Dii l°L\ + D h D h (-^-) - - (1.2) 
og a b yogab^J \Ogab,hi 2 ) 

of some locally defined Lagrangian 

L — L(x , g c d, gcd,ii j gcd,iii 2 i ■ ■ ■ ) gcd,i\i 2 ---i x ) 

depending on the components of the metric tensor and their derivatives. Here Di 
denotes the standard coordinate total derivative operator. 
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If the Lagrangian L transforms as a scalar density, that is, 

£ prX5 L = prX 5 (L) + div£L = 0, for all X ( e g, (1.3) 

or, equivalently, the density L g = g^^L is invariant under the prolonged action of g, 
then the Euler-Lagrange expressions T ab = E ab (L) constitute a tensor density, whereby 

£ P rxJ ab = pr X^(T ab ) - 2&T b)c + dw£T ab = 0, for all X 5 e g. (1.4) 

Here C denotes the standard Lie derivative operator. Lagrangians L and differential 
operators T ab satisfying (11. 3p or (II. 4p are also known as natural tensor densities. 

In addition, in light of the diffeomorphism invariance (ll.4p . Noether's second theo- 
rem implies that the components T ab = E a6 (L) are divergence-free, 

D b T ab + r£ c T 6c = 0, (1.5) 

where the T^ c denote the standard Christoffel symbols of the metric g. 

In the present paper we investigate a partial converse to Noether's second theorem 
for third order operators for metrics, that is, whether a symmetric, type (2,0) differen- 
tial operator T ab = T (ab \x\ g cd , g cdM , g cdjili2 , g cdAli2i3 ) satisfying the invariance condition 
(11.41) and subject to the differential constraints (I1.5P coincides with the Euler-Lagrange 
expression of some Lagrangian L. 

Our main result is the following. 

Theorem 1. Suppose that a third order differential operator 

rriafo _ rp(afr)/ i \ 7 _ -| 

1 — 1 yX , g cd , g c d,ii i gcd,iii2 i gcd,iit2i3 ) ; fl, — 1, . . . , 771, 

admits the symmetries (11.41) and satisfies the differential constraints (11.51) . Then T ab 
is locally variational. Moreover, suppose that T ab is everywhere smooth. Then it can be 
written as 

Tab _(E ab (L), ifm = 0,1,2 mod4 ; 

\E ab (L) + C ab } ifm = 3 mod 4, 1 ' ' 

where the Lagrangian L is a scalar density satisfying (II. 3p and Cp b is the generalized 
Cotton tensor associated with an invariant polynomial P on So(j>, q), p + q = m, of degree 
(m+l)/2. 

The generalized Cotton tensors Cp b are, therefore, locally variational, but, as is 
proved in [5], they can not be written as the Euler-Lagrange expressions of a natural 
Lagrangian. In the physically most relevant situation with m = 4, Theorem [1] asserts 
that if a third order, natural system of differential equations 

T ab (gcd, g c d,h , 9cd,nn > 9cd,nni 3 ) = °, where T ab = T (a6) , 
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for the components g c d of the metric tensor is divergence free ( 11. 51) . then there is a 
natural Lagrangian L so that T ab = E ab (L). 

This paper is organized as follows. After covering some preliminary material relevant 
to the problem at hand in section [2j we analyze in section [3] the relationship between 
symmetries (11.4j) and the conservation law (11. 5p for metric field equations. In particular, 
we show that any natural differential operator T ab admitting translational conservation 
laws is necessarily divergence free. This interesting though elementary fact does not 
seem to have been previously noted in the literature. Then in section H] we present 
the proof of Theorem [1] and, finally, in section [5] we discuss some open problems and 
generalizations of the work at hand. 



2 Preliminaries 

In this section we collect together some basic definitions and results from the formal 
calculus of variations on jet spaces germane to the problem at hand. For more details 
and proofs we refer, e.g., to [2j [23] . 

Let G -> M m be the trivial bundle of metrics, that is, of non-degenerate sym- 
metric bilinear forms on IR m with fixed signature. Denote the coordinates of R m by 
(x 1 , x 2 , . . . , x m ). Then the components g ab = g^ ab ) of a metric g are determined by 
g = g a bdx a <8> dx b , so that, as a coordinate bundle, 

G = {(x\g ab )} -> {(V)}, where a < b. 

We denote the bundle of order k jets, < k < oo, of local sections of G by J fc (G); in the 
induced coordinates 

J (G) — { (-^ ) 9ab) 9ab,i\ i 9ab,i\%i i ■ ■ ■ j 9ab,i\%i---%i > • • • i 9ab,i±i2---ik )} ' (^"^) 

where g a b,iih-ip a < b, stands for the l-th order derivative variables. For notational 
convenience we let g a b,hi 2 -ii = 9ba,hi 2 -iii when a > b. We also use to collectively 
designate all the variables g a b,i r -i p , V - 0, . . . , k, up to order k. 

Let / = (^1,22, • • ■ 1 < i p < m, denote an unordered multi-index of length |/| = /. 
Define partial derivative operators <9 a6,/ , 1 < a, b < m, \I\> 0, by 

fx axb ^s i2 —8 ik) if l/l = I 71 

& xb ' I g* i j= (c d) jl h 3k ' 11 1 (2.2) 

where round brackets indicate symmetrization in the enclosed indices. Then, for exam- 
ple, the standard coordinate total derivative operators Di on J°°(G) are given in terms 
of the differential operators (12.21) by 

A = ^7 + Y,9ab,nd abJ , * = l,2,...,m. (2.3) 

|J|>0 
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The expression ( 12 .31) leads to the commutation formula 

[& ib ' 1 , Dj] = " h df } . (2.4) 

We will employ the standard Einstein summation convention in what follows. 
The flow of a vector field 

X = P\ x \ 9cd )^- + Q ab (x\g cd )d ab (2.5) 

on G induces a transformation on the space of sections of G, and, hence, by differen- 
tiation, it generates a local 1-parameter transformation group acting on J fc (G), k > 0. 
The associated infinitesimal generator is called the k-th order prolongation of X and 
is denoted by pr k X. The components of pi k X are given by the usual prolongation 
formula 

W k X = P l D t + £ D I (X w , ab )&*> 1 , (2.6) 

|7|<fc 

where the X eVjab denote the components of the evolutionary form 

X ev = (Qab ~ P C 9ab,c)d ab 

of X and where, given a multi-index I = . . . we use the abbreviated notation 
Di = Di 1 ---Di k . We will also write pr°° X = prX. The vector field (12. 5p is called 
projectable if the coefficients P a = P a (x l ) are functions of the independent variables x 1 
only. In particular, the infinitesimal generators of the action by the lifted diffeomor- 
phism group (II. ip form a Lie algebra g of projectable vector fields on G with 

Xt,ev,ab = ~^i%9b)c ~~ C9ab,c- 

We associate to a given differential operator T ab = T ab (x\gW) the source form 

T = T ab dg ab Av, (2.7) 

where v = dx 1 a dx 2 a • • • a dx m is the standard volume form on M. m . The source form 
T is called natural if T is invariant under the prolonged action of the diffeomorphism 
group g, that is, 

£ prX? T = 0, for allf e X(R m ), 

where X^ is defined in (11. ip . As is easy to verify, a source form T = T° dg ab hv is natural 
precisely when the components T ab form a tensor density as in (jl.4p . A vector field X 
on G generates a conservation law for T if there are differential functions t l = t l (x^ , g^), 
i = 1, . . . ,m, so that 

A" ev ,abT ab = Djt\ (2-8) 

The source form T is said to be derivable from a variational principle if there is a 
Lagrangian function L = L(x t ,g^) such that T is the Euler- Lagrange expression of L, 
i.e., 

T ab = E ab( L j = £ (_ jD ) / (^6,/ L ^ ( 2 9 ) 

|/|>0 
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It will be convenient to call 

A = L(x a ,g [l] )u 

a Lagrangian n-form and 

E(A) = E ab (L) dg ab a v 

the Euler- Lagrange form associated with A. As is well known [2], the Euler-Lagrange 
operator commutes with the prolonged action of projectable transformations on G; 
infinitesimally, 

E(£ prX A) = £ prX E(A), (2.10) 

for every projectable vector field X and Lagrangian form A, where C denotes the stan- 
dard Lie derivative operator. 

The Helmholtz operator Ht acts on evolutionary vector fields Y = Y a f,d ab on G by 

H T (r) = £ pry T-E(r- T), (2.11) 

cf . [9] . If we write 

H T (y) = Y, (DjY^R^dga, a u, (2.12) 
|/|>o 

then the components H T b ' cd,/ of Ht are explicitly given by 

R af>,cd,I = d cd,I T ab _ (_ 1 )|/| E ab > J( T c^ ) |j| > q 

Here the E a6,/ denote the higher Euler-Lagrange operators [2] acting on a differential 
function F defined on some J k (G) by 

E ab,i (F) = £ (WM)(-D)j(ff*» IJ F), \I\ > 0. 

|,/|>0 

Note that if T = T ab (x l , g^)dg a b a v is of order k, then H T b ' cd ' 7 = for |/| > k and for 
\I\ = 0, . . . ,k, the components H T 6,cd ' 7 are of order at most 2k - |/|. 

It is not difficult to see that a source form T = E(L) deriving from a variational 
principle satisfies the Helmholtz conditions Ht = 0, or, in components, H T 6 ' cd,/ = 0. 
Conversely, one can show [2] that if the Helmholtz conditions Ht = are satisfied, then, 
at least locally, that is, in some neighborhood of each point in its domain in J fc (G), 
the source form T can be written as the Euler-Lagrange expression of some Lagrangian 
L. Accordingly, we will call a source form satisfying the Helmholtz conditions locally 
variational. 

Proposition 2. Suppose that X = P t d/dx l + Q a bd ab is a projectable vector field on 
G and that a source form T is invariant under the prolongation prX of X . Then the 
components H T 6 ' cd,/ of the Helmholtz operator 'Hp associated with T satisfy the invariance 
conditions 

r)P3 

pr Xm^' 1 ) + V m^cT^Qefj + Yd' cdJ d ab Q ef + ^—YL^ 1 = 0, (2.14) 

where Q c d,j denotes the g^j -component ofprX. 
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Proof. We first compute 

r )) = r v/7'. > - f... v WY _. 

(2.15) 



£ prX (H T (y)) = £ prX (£ pry T) - £ pr xE(F - T) 



= £ pr[ x,y]T - E((£ pvX Y) - T) = H T ([X, F]), 

where we used the invariance of T and the equivariance of the Euler-Lagrange operator 
under the prolonged action of projectable transformations. Next write Y = Y cd d cd . 
Then, on account of (12.121) . 



£ wX (K T (Y))= £ ( w X{B.f cd > I )D I Y cd + Kf cd > I w X{D I Y cd ) 
|/|>o V 



+ Rf^'iDjY^Qef + ^- R^' 1 DjY cd \d gab a v. 



(2.16) 



The identity 

[prX,prF] =pr[X,F] 



yields 



VTXiDjY*) = Dj([X,Y] cd ) +V iY(Q cdJ ) 

= Dj([X, Y] cd ) + Y, (DjY ef )d^ J Q cdJ . (2-17) 

|J|>0 



Now by virtue of (I2.17p . equation (I2.16P becomes 

£ piX (R T (Y)) = £ (vrX(E^)D I Y ei + B^D I ([X,Y] ei ) 



+ ^ ^\DjY ef )d^ J Q cdJ + Kf^iDjY^Qef (2 18) 

[J|>0 



9P 3 u ab,cd,I 
~Z — r ilrp 

dx3 T 



H°' c ' DjYd \dg ab a z/. 



A comparison of (12.151) with (12.181) yields equation (I2.14p . as required. \qed\ 



ab,cd, 
T > 



Proposition 3. Let T be a third order source form. Then the components H 
H^' ,l ; K^,' ,iJ , FT^' ,iJ of the Helmholtz operator associated with T satisfy the inte- 
grability conditions 

YLf cd + R^ ab = DiE^' cd4 - /;,/), I r;''-"'- ; ' + D i D j D k l£*' cd ' ijk , (2.19a) 
JJ06.0M _ R Mi = 2 D j R^' cd ' ij _ 3DjDk}1 a T b ^\ (2.19b) 

R aJ>,cd,ij + jjaJ.oft.ii = 3jDfeH ^' Cd ' ijfc , (2.19c) 

H« 6 - cd « fe -H^' aMife = 0. (2.19d) 
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Proof. Equations ( 12.191) can be derived by direct computation using the expressions 
(I2.13P for H T 6 ' c<i ' 7 . Since T of order 3, these yield 

-^ab,cd,ijk _ Qcd,ijkr^ab _|_ Qab,ijk<-£cd ^2 20) 

which immediately implies equation ( I2.19d[) . Likewise, by ( 12.13p . 

E^' cd ' ij + H^ 6 ^' = 3D k (d ab > ijk T cd + d cd ' ijk T ab ) = 3D k R^ b ' cd ' ijk , 
so that ( I2.19c[) holds true. Next we compute 

Hab,cd,i -rjcd,ab,i 
rp — rlrp 

= 2Dj(d cd ' ij T ab - d ab > ij T cd ) + 3D j D k (d ab ' ijh T cd - d cd ^ k T ab ) 

= 2D j (d cd ' ij T ah - d ab > ij T cd + 3D k d ab ' ijk i: cd ) - 3D j D k {d ab > ijk T cd + d cd > ijk T ab ) 

= 2D i H£' c<Mj - 3D j D k R a ^ cd ' ljk , 

which yields (12. 19b[) . Finally, 

H T b ' cd + R cd ' ab = A(<9 C(M T a6 + d ab ' l T cd ) 

- D l D j (d cd ' ij T ab + d ab ' ij T cd ) + D i D j D k (d cd ' ijk T ab + d ab ^ k T cd ) 
= D i (d cd ' i T ab + d ab '*T cd - 2 J D j a ab ' ii T cd + 3D j D k & ih ' i: > k T cd ) 

- D l D j (d cd > ijr T ab - d ab > ijr T cd + 3D k d ab > ijk T cd ) + D i D j D k {d cd > ijkr T ab + d ab ^ k T cd ) 
= M ' 1 - DiDjEf ^ + DiDjD k R^ ,cd,ijk , 

which proves the first identity (I2.19ap . \qed\ 

Remark 4. The Helmholtz operator Ht can also be characterized in terms of the 
differential 5y of the variational bicomplex [2] as Ht = <5yT, whereby the foregoing 
conditions ensue from the general identity 5y = 0. 

The following Lie derivative formula, as established in j2j[7], is central in the proof 
of our main Theorem. 

Proposition 5. Let T be a source form and X a projectable vector field on G. Then 

£ piX T = E(X CV -i T) + H T (X ev ). (2.21) 

An extension of the Lie derivative formula ( 12.21 j) to non-projectable, generalized 
vector fields can be found in [2] . If T is a locally variational source form, then equation 
(I2T2TD reduces to 

£ P rxT = E(X ev -i T). 

Now 

£prxT = 0, 
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if the source form is invariant under X, that is, if X is a distinguished symmetry of the 
system T ab = 0, while 

E(X CV ^T) =0, 

provided that X generates a conservation law for T; see [2]. Thus equation (I2.2ip fur- 
nishes a version of the classical Noether's theorem for projectable vector fields expressed 
directly in terms of the system of differential equations without explicit reference to a 
Lagrangian. 

On the other hand, in the situation of Takens' problem, each vector field X belonging 
to the Lie algebra T is, by prescription, a symmetry of the source form T and generates 
a conservation law for T. These requirements lead to the conditions 

H T (X ev ) = 0, for all X e T, 

on the Helmholtz operator of T. The primary objective in the analysis of Takens' 
problem is to identify, on mathematical or physical grounds, interesting classes T of 
source forms (i.e., differential equations) and symmetry algebras T of vector fields so 
that one will be able to classify all T-invariant Helmholtz operators Hx with T e T 
satisfying the conditions Ht(X 6v ) = 0, X e T. 

We conclude this section by briefly recalling the construction of the generalized 
Cotton tensors C p . For more details and proofs we refer to [5]. 

Let u)j = r i j k dx k , i,j = l,...,m, denote the connection forms of the Riemannian 
connection associated with a metric g, and, as usual, write 

fij = dUj + Ulj. A Uj 

for the associated curvature form. In terms of the components of the Riemannian 
curvature tensor, 

= -R i jkl dx k a dx l . 

Next let 

tr r :gl(m) -> R, tr r (A) = trace(A r ), r = l,2,..., 

denote the elementary polynomials on the Lie algebra gl(m) invariant under the adjoint 
action of the general linear group Gl(m). As is well known [22], any 0(p,q), m = p + q, 
invariant polynomial P € I 2n (so(p,q)) of homogeneous degree 2n on the Lie algebra 
So(p,q) of the orthogonal group 0(p,q) can be expressed as P = p(tr 2 ,tr 4 , . . . ,tr 2s ) for 
a uniquely determined polynomial p:M s -> R. 

Next assume that m = An- 1. We associate to a polynomial P = p(£r 2 , £r 4 , . . . , tr 2s ) 6 
J 2n (so(p, q)) the second order differential functions Ep bc (g) by 

E£ bc (g)v = g ad P c d (n) Adx\ 

where the P^ are the components of the derivative of P, that is, 

lp(A + tB)\ t = = P c d (A)B d c . 
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Now the generalized Cotton tensors are defined by 

C P (g) = C?d gab Au, 

where the components Cp b are, by definition, given by the covariant derivatives 

Cf = ^ c (E^ cb (g) + E b P ca (g)). 

One can show that C p (g) is a third order, locally variational natural source form [5J. 
When m = 3 and P = tr 2 , the Cotton tensor C p agrees with the well-known Cotton- York 
tensor, [121 ED]. 

The following theorem is proved in [5j. 

Theorem 6. Let T be a natural, locally variational source form. Then T can be ex- 
pressed as 

fE(A), ifm = 0,l,2 mod 4; 

■E(\)+C p (g), ifm = 3 mod 4, 

where A = Lv is a natural Lagrangian form and C p is the generalized Cotton tensor 
associated with a uniquely determined Pe I~2~(so(p,q)). 



3 Symmetries and conservation laws 

In this section we analyze the relationship between the diffeomorphisms symmetries 
(jl.4p and the differential constraints (jl.5p expressing divergence-freeness in metric field 
theories. As is well known, the commutation formula (l2.1Up . together with Noether's 
theorems, implies that the Euler-Lagrange expression T = E(A) of a Lagrangian form 
A with symmetries (11. 3p possesses symmetries (11. 4p and is, in addition, constrained by 
the identities (II. 5p . The following result, which bypasses the Lagrangian and is stated 
directly in terms of the system of differential equations, is a slight but non-vacuous 
extension of the above conclusions furnished by the classical Noether's theorems; see 

Proposition 7. Suppose that the source form T = T ab (g^)dg a bAi/ is locally variational. 
Then T is natural, £ pr x e T = for all e q, if and only if it is divergence- free, 
£» 6 T a6 + r» c T 6c = 0. 

Proof. By assumption the Helmholtz operator Ht of T vanishes, and so equation (12.211) 
reduces to 

C prX J = E(X € , ev T) = -E hk (2^ a g bc T ab + ^9^)6^ a v. (3.1) 
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Recall that the Euler-Lagrange operator annihilates total derivative expressions, E(D a F) = 
for all differential functions F = F(x l ,g^) and a = 1, ... ,m. This permits us to inte- 
grate the right-hand side of ( 13.1 j) by parts to conclude that 

V h \2^ a g b X b H c g ab ,X b )dg hk Av 

= -E hk (^(2g bc D a T ab + 2g bc , a T ab - ^ c g ab , c i: ab ))dg hk a v 
= -2E hk (^g bc (D a T ab + T b cd T cd ))dg hk a v. 

Hence 

C piX J = 2E hk (£ c g bc (D a T ab + T b cd T cd ))dg hk a (3.2) 

which immediately establishes the Proposition in one direction. 

It then remains to prove that the condition £ pr x e T = for all £ € X(M. m ) implies that 
D a T ab +T b cd T cd = 0. For this, suppose that for some index b a , D a T ab ° +r^T cd = F(x\ gW) 
is of order I and that for some h, k and J with | J| = i, 

{&*> J F)(4,gP)*0, (xlgi l] )e.J l (G). 
Now choose £; e ^(R m ) such that 

-Q^f(x l o)9o,bc = 5 bbo , -^f(x l o )g oM = 0, K + J. 

Let (xj,,#P ] ) e J 2 '(G) be any point projecting to g\ . Then 
E /lfe (^ c (^aT ab + r^T cd ))(^,^ ] ) 

= (-1)11^ (^(^(^.t- 6 + r^T cd ))) ] ) = (^ fe - / F)(4,^ ] ) * o, 
\I\<1 

which is a contradiction. Thus D a T ab + T b cd T cd = h^x*) are functions of x % only. But 
due to translational invariance of the source form T, each h b must be constant. Finally, 
by the definition of the total derivative operators (12. 3p and the Chrisoffel symbols T b cdl 
the expressions D a T ab + T b d T cd vanish when evaluated on the jet of the constant metric 
g = diag(l, . . . , 1, -1, . . . , -1), showing that h b = 0. 511 

Recall that a source form T = T ab dg ab a v admits translational conservation laws if 

E(g ab)P T ab ) = 0, for all 1 < p < m. 

Proposition 8. Suppose that a source form T is natural and admits translational con- 
servation laws. Then the covariant divergence o/T vanishes, D b T ab + T?T bc = 0. 
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Proof. We will prove this Proposition by showing that 

H T (JQ jCV ) = for all f e X(R m ). (3.3) 
First, with X = X ( = £ a d/dx a - 2^ a g bc d ab e 0, equation (1213)1 becomes 

prX ? H? cd ' 7 + J) (a cd ' / X^ e/)J )H^ e/ ' J - 2e, ( e a HT )e,cd,/ + ^ e e H^ ),cd ' / = 0, (3.4) 

where ^, e /,j = P r -^Q(^e/,j)- The source form T is also invariant under r p = d/dx p , so 
the hypothesis and equation (I2.2ip imply 

£ gcdjp&t *' 1 = 0, for all 1 < p < m. (3.5) 

|/|>o 

Now apply the vector field prX^ to the above equation to see that 

£ .Y i .,.,, /; ,li:;'"" u + £ ^, /p prX € (<^ 7 ) = 0. (3.6) 
|/|>o |/|>o 

Combining (13.4jl and (13.61) we obtain 

|/|>o |/[>o \[J|>|/| , 

E v -rjab,cd,I r> ^ / Qcd,I y \Tjab,ef,J _ n 

^5,cd,7" P -tl T - 2j L 9cd,l P \0 A^ >e /,jJil T - U, 

|/|>0 |/|>0|J|>|/| 

where we used (13.51) twice. On account of the definition of the total derivative operators 
(12.31) . the above equation can be written as 

r) X 

E v jjab.cdj n v Tjab,cd,I uy *-£,cd,I jjab,cd,I _ n to 1\ 

|/|>o |/|>o |/|>o ax 

By the standard prolongation formula (I2.6p . 

X^cdjp = D p Xt: >cd j - —gcdjq- 
Thus equation (13. 7p simplifies to 

Ed? jjab,cd,I , V 1 ®X£ t cd,I jjab,cd,I 
-^~:9cd,iqn- T + tt~: n T 

\I\>0 OX \I\>0 OX 

oi; m t - Z, ^^/axP, C d,/W T - u, 

|/|>0 ux |/|>0 
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where we again used ( 13.5H . The vector field £ e X(M. m ) is arbitrary, which allows us to 
conclude that 

£ X^jH^ 1 = 0, for all f € <Y(R m ). 

|/|>o 

But on account (13.51) . 

E v Trab,cd,I _ / n v cq \jrab,cd,I _ r> n y Traft,cd,/ _ n 

^£,cd,ltL T - 2^ K-LJl^^ev^d + C, 9cd,Iq)n- T - ^ JJl^£,ev,cd£i T - U, 

|/|>0 |/[>0 |/|>0 



that is, equation (13. 3 j) holds. 

Due to the g-invariance of the source form T and condition (13.31) . the Lie derivative 
formula (I2.2ip now yields 



m c 9bc{D a T ab + T b cd T cd )) = 0, for all £ c € X(R m ). 

Next we continue as in the second part of the proof of Proposition [7] to conclude that 
the source form T is divergence- free. \qed\ 

4 Proof of Theorem 1 

The proof of Theorem [1] relies on the following Lemma, which is a special case of a 
more general result presented in [61 [22] . 



Lemma 9. Let T = T ab dg ab Ais, T ab = T (ab) (x\ gW), be a k-th order source form, where 

a r 
bc- 



k > 1, and assume that the covariant divergence D b T ab + T? T bc = vanishes identically 



Then the component functions T ab are polynomials in the k-th order derivative variables 
9cd,h-ik °f degree at most m-1, where m is the number of the independent variables. 

Proof. By assumption, 



dT ab 



dx b |,| n 



+ X 9cd,ibd cd ' Ir T + T a cd T cd = 0. (4.1) 



Now terms in (14. ip involving the order k+1 variables g a b,j, \J\ = k + 1, yield the equations 

d cd,(I T b)a = Q ( 4 _ 2 ) 

Write c^f x for the partial differential operator 

dfx = E Xi^' 1 = Z dah ' il " 4kX ii-Xi k , (4.3) 
\i\=k 

where X = {Xi, . . . , X m ) e (R m )* is a covector on R m . Then equation (14.21) is equivalent 
to 

X b d c k d x T ab = for all X. (4.4) 
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Next let G a ' Cldl "' Cmdm denote the mappings 

Qa,cidi-- c m d m ^ j£l X m Y") = Q Cldl . . .Q c "idm'-£abY^ 

The operator T ab is polynomial in g a b,i, \I\ = k, of degree at most m - 1 if and only if 
all the mappings G a ' Cldl "' Cmdm vanish identically. But by ( 14.4ft and linearity in Y, the 
equation 

G a ' Cldl "- Cmdm (X 1 ,...,X m ,Y) = 

holds whenever F is a linear combination of the covectors X 1 , . . . , X m . Consequently 
Q^adv-cmdm van i s hes for almost all X 1 , . . . , X m , y e (R™)*. By continuity, G f0 > Cldl "" Cmdm 
must vanish identically. [all] 

We next employ the Lie derivative formula ( 12.21)) to derive key identities among the 
components of the Helmholtz operator associated with a third order natural, divergence 
free source form. On account of these identities and the integrability conditions of 
Proposition [31 the proof of the first part of Theorem [1] reduces to showing only that the 
third order components of the Helmholtz operator vanish. 

Proposition 10. Suppose that T = T ab dg a b a v is a third order natural, divergence free 
source form, that is, 

C piX J = 0, for all X € eg, and D b T ab + T^ c T bc = 0. 

mi j.7 j. TT(ib,cd TT(ib,cd,i -nab,cd,ij T-rab,cd/ijk r ±i tt i i i± ± tt 

1 hen the components H T , H T , H T , H T of the Helmholtz operator Ht 

associated with T satisfy the following identities. 

Tjab,cd TTab,cd,i Trab,cd,ij Trab,cd,ijk _ n /» r \ 

+2< ?cd , ie H? C(i ' i/ + 2g cem R ab ^ k + 3g cd>ije R$' cd > ijf = 0, (4.5b) 
2g ce K a ^ + 4g ceJ Rf cid > i)j + g cf , e K$' cf ' di 

+ 6g ce , k Rf c{d ^ k + SgcjjS^^ = 0, (4.5c) 

2< ?ce H? c( ^ ) + 6fe, fc H? c(( ^ )fc + ^ ie H^' cfc '^ d = 0, (4.5d) 

Proof. By the Lie derivative formula (12.211) . 

H T (X c , ev ) = (DjX^ cd )R^ cdJ dg ab A v = 0, (4.6) 

for all € g, where 

X ?>cv = X^ cv , cd d cd = -(2^ (c g d)e + i e g cd ^d cd . 
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Let /i cd ' 7 = yu( cd )'( 7 ), |/| > 0, be multivectors. We compute 

I I t| \^cl9de,jl-' CC 

\il\J\>o y l J l ; 

= EC/E( |J| |t^r%^ 1 '^ 

|/|>0 [J|>0 v l J l 1 7 



^(Ccfe)^' 7 = E ( |J| l + J J| )Cc/^,^ c 

|/LUI>o v l J l 7 



(4.7) 



and 



^(e^, e )^ 7 = E E (' ;J 'k^ cd,/J - ( 4 - 8 ) 

|/|>0 |J|>0 v l J l ' 

The derivatives £ e 7 are independent, so in light of (14. 7p . fl4.8j) . the coefficient of £ e 
in (14. 6p yields the equation 

E^e<' cd ' J = 0, (4.9) 

[J|>0 

while for III > 0, we obtain 



E 



( |/| l ; l ' l ;> fe ,,Hr (fa *-- , ' ,j + ( |/| l + / l )^ J «H 



ab,cd,IJ 



\j\>al' 



0. (4.10) 



Keeping in mind that H T b ' c<i ' J = for 1 1\ > 4, equations (14. 5 p follow from (14. 9 p and (I4.10p 
by inspection. \qed\ 

We note that if n abdc = ^( ab )( cd ) is a valence 4 tensor satisfying the cyclic identity 
^ a (bcd) _ 0> then 

3 

^abca! _ ^cdafe _ _ ^^a(bcd) _j_ ^b(acd) _ ^c(abd) _ ^d(abc)^ _ g ^ -^-Q 

so that ^ abcd is symmetric under the interchange of the pairs of indices ab and cd 

Corollary 11. Suppose that T = T ab dg a f, a u is a third order natural, divergence free 
source form. Then T is locally variational provided that the third order components 
jja6,cd,y Q j ^ g Helmholtz operator Ht associated with T vanish. 

Proof. We shall show that as a consequence of the assumptions, all the components of 
the Helmholtz operator Ht vanish identically. The proof is based on identities (I2.19P 
and (14.51) of Propositions [3] and [101 respectively. 

We first observe that by (12. 19cj) . the second order components H™' are skew- 
symmetric under the interchange of the index pairs ab and cd. Due to the assumptions, 
equation (I4.5d[) reduces to Y^' c ^ d ' 1 ^ = 0. Consequently, identity (14. lip applied to the 
index pairs cd and ij shows that the Yi c ^' cd,l: ' are also symmetric under the interchange 
of these pairs. Thus the second order components H T 6,cd ' 1 "' must vanish identically. 

By the above, equation (I4.5cl) now implies that fj^' c w) = o. Thus }^' cd ' 1 i s symmet- 
ric and skew-symmetric in overlapping pairs of indices and hence must vanish. Finally, 
equation (14.5b|) presently shows that H T b,cd = for all a, b, c, d. \qed\ 
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We will use a semi-colon to indicate differentiation with respect to the weighted 
partial derivative operators ( 12 .2)) . so that, for example, 

Qab,ij Qcd,klm jp _ p\ab,ij\cd,klm 

for a differential function F on J°°(G). 

We will use the following result in the course of the proof of Theorem HJ 

Proposition 12. Let T = T ab dg ab av be a third order divergence-free source form and 
let h^ 6 '^'* 1 * 2 denote the second order components of the Helmholtz operator associated 
with T. Then 

(i) r Y ab \ c i c 2,iii2(ki\;did 2 ,\jij2j3) _ r£i 1 i 2 ;c 1 C2,ab(i 3 \;d 1 d2,\jij2j[i) . f£ 

(il) r-p(a6|;ciC2,|ii)i2'j3 _ ^fj«2«3;ciC2,a6ii . ^ 

3 

(m) ^ ab ,ciC2,iii2;did2,jij2j3j4 _ ^<j,iii2;ab,ciC2(ji\;d 1 d2,\j2j3j4) ^ 

where indices enclosed within vertical bars are omitted in the symmetrization. 
Proof. Since T is divergence free and of order 3, equation (14. ip yields 

rjia(b|;cic 2 iiii2«3) _ Q (A -^tj\ 

Consequently, 

rjna(fe[;CiC2|jl«2)'3 _ _ _ fja«3 ;ci c 2 ,bil %2 (A ~^Q^ 

3 

which, on account of (I4.15p . implies that 

rpa(6|;c 1 C2|iii2)(«3|;did2,b'ij2i3) _ _^^a(iz\\c\C2&\i-2\>i\d.2\d\hh) 

3 

_ _l^CiC 2 ,feuj2 r J i a(*3l;dlrf2,|ilj2i3) _ Q 

3 

Equation (fl~TTj) now yields ( l4~T2j) . 
Next, due to (14.151) . we have that 

rj(a*[;ciC2,|ii)i2i3 _ _gy m r J 1 (*2|a;ciC2,feii|«3) _ ^ rj«2i3;ciC2,a&ii ^ -^ijA 

{afeii } 3 

where in the second step we used equation (14.161) . 
Finally, in light of the commutation formula ( I2.4p . 

-^ab,ciC2,iii 2 ;d 1 d 2 ,jij2j3j4 _ Qd\d2 ,jij 2 jaj4 (QciC2,hi2<~£0- b _ ^ab,n«2 r J 1 ciC2 +Q£). ^^a6,iii2«3 r p c i c 2 ^\ 

= 3 r p clC2;a ^'* 1 * 2 ^ 1 ' ;dl ^ 2 '- ?2 - ?3 -' 4 ^ = 3 r Y iii2 '> ab ' CiC2 (j i \'> did2 \i 2 33j4) 
which completes the proof of the Proposition. \qed\ 
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Proof of Theorem d Expression (II .61) is a straightforward consequence of Theorem [6] 
for locally variational source forms. Thus it suffices, by Corollary [TTJ, to show that the 
third order components H^' ' y of the Helmholtz operator Ht vanish. 

In order to streamline our notation, we let p = (p±,P2) and q = (qi,q2, 93) denote 
multi-indices of integers 1 < p t < m, 1 < q t < m of respective lengths 2 and 3, and write 
$p,q = Qv\V2A\i2Vi j n w hat follows. We also employ the notation 

^;piqi;p 2 q2; -;p Is q B _ ^pi ,qi ^P2,q2 . . .^p s ,q s p 

for the repeated derivatives of F with respect to the weighted partial derivative opera- 
tors d Pu,ciM = 5 PulPu2 ' 9 " 19 " 29 " 3 . 

By Lemma [9], the m-fold derivatives T ab,Piqi ' P2q2 ' "' PmC|m = 0, and consequently, the 
third order components Y['^> cd ' l i k = ^cd.jjfc^afr^aZw'fcx^ G f ^g Helmholtz operator satisfy 

^ab,cd, ijfc;piqi;p2q2;--;Pm-iqm-i _ g 

We will prove by induction that H^* i * ;piqi:paqa; "" !Pp,,r = for aU < r < m - 1. 
In order to carry out the induction step we will assume that 

H a6 J cd,iifc;p 1 q 1 ;p 2 q 2 ,.. i p r+1 q r+1 = ^ for some Q < r < m - 2. (4.18) 

Our first goal is to show that 

pja6,cd;Zii2,nin2n3n4;piqi;P2q2;— JPrqr _ g ^ ]^g a ^ 

jjo6,cd,i;ZiZ2,nin2n3n4;piqiiP2q2;- , -;Prqr _ g ^ 19b) 

^ab,cd,ij;l 1 l 2 ,n 1 n 2 n 3 n4;p 1 ci 1 ;p2(l2; --,PrClr _ g ^ ^g c ^ 

For this, we start by applying the operator d lll2 > nin2H3n4 to (12. 19cj) to see that 
Thus by the induction assumption (I4.18p . 

■^ab,cd,ij;l 1 l 2 ,n 1 n2n 3 n4-,p 1 ci 1 ;p2<l2-,- --,PrCir + j^cd,a6,ij;ZiZ2,ni»i2"-3™4;Piqi;P2q2;- -;Prqr _ g 

that is, the components }^' cd > l i' lil2 > n ^ n2n 't ni are skew-symmetric under the interchange 
of the index pairs ah and cd. On the other hand, repeated differentiation of (14.5dl) 
yields the equation 

jja6,c(d,ij);ZiZ2,nin2ri3n4;piqi;P2q2i"-;Prqr _ g 

Thus by (14.111) . the components Y^' cd)l: '' lll2,nin2nAn " 1 are symmetric in the pairs cd and 
if. It now follows that 1 )4.19 eft holds. 
By differentiating ( I4.5c[) we see that 

jjofe,c(ci,i);ZiZ2,nin2".3n4;piqi;P2q2;- -;Prqr _ g 
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where we used (I4.19cl) and the identity Y{^ cd ^ k ' l ^' n ^ n ^ n i = q Consequently, ( I4.19bl) 
holds. Finally, the first equation (14.1 9ap now follows from (I4.5bj) . again by differentia- 
tion. 

Next apply the differential operators d lll2 ' nin2n3n4 , <9 Pl,qi ,. . . ,5 Pr,qr , to equation (|4.5aj) . 
Keeping in mind that the components Y{^' cd ^ k are of order 3 we conclude that 

TTab,cd;W2,n.in2n3n4;piqi;---;Prqr , TTa6,cd,i;W2,niri2n3n4;piqi;-";Prqr 
ycd,e rl - r £ ' 'Jcd,ie rl - r £ 

, TTo6,cd,tj';hia,mfi2n3ri4;piqi;*";Prqr 



(>H\ ■[]"''-;': i,»/j-|<// »//:: /■.../(!"-".;" i :pjQ J :-p/q, ; -:p, q, (4.20) 

t=l 

+ ^(ii|jja6,iiZ 2 ,|n2n3n4);piqi;— ;Pr-q r _ g 

where the hat indicates omission. On account of (14. 191) . the above equation reduces to 

r 

h 6 T (4.21) 

_l_ ^(ni|jjo6,2ii2,]n2n3n4);piqii---;Prqr _ g 

Our next goal is to show that H^' cd ' e ^'' 1 ' 2 '™ in2 ™ 3n4 ' Piqi, '''' Pr ~ iqr ~ 1 . As a result, equation 
( 14.21 j) becomes 

^(ni|jja6,ZiZ2,|n2n3n4);piqi;— ;Prqr _ g ^ 22) 

which will immediately imply the induction step. 

To proceed, we first contract in the index pair e, qu in (I4.2ip . which yields the 
equation 

m + _ jj-o6,piiPi2,gi2gi3^ib,nin 2 n3n4;p2q2;"-;Prqr 

3 T 

r 

+ ^ JJ a6 'PtlPt2,(<?tigi2|;hi2,nin2n3™4;Pl,|9t3)9l2<7l3;P2q2; --Ptqi --;Prqr (4.23) 
t=2 

_l_ jjafe,/i« 2 ,(nin2n 3 |;pi,|n4)(jri2(3i3;P2q2;---;Prqr _ g 

On account of (14.141) . we have 

jjafe,pnpi2, gi29i3 ;'i'2,ni«2"3"-4;P2q2;---;Prqr _ gfp?i29i3;<i6iPiiPi2(ni|;/ii2,|n2n3n4);p2q2;- , -;Prqr 

We next use (I4.13P and (14.141) to compute 

jja6,ptipt2,(9ti9t2|;iii2,nin 2 Ti3n4;pi,|9t3)gi29i3;p2q2;— ptqt— ;P7-qr 

= 3 Sym T^ tl9 * 2 ' ;ab ' Ptip * 2 ™ i; ' 1 ' 2 '™ 2n3n4;pi ''' 3t3 -'' 3l2 ' ?13;P2q2; " 'P^q*""P'-q'- 
{11-14} 

= 3T(* 1 ' t2 l'P 1 'l* 3 )5 12 9 13 ' a ^P' 1 P' 2 (" 1 li' li2 'l n2n3n4 )'P 2 ^ 2 "'P^'''!P r( l'' 

_ rpgi2gi3;piiqt;a6,ptiPt2(ni|;(iZ2,|n2n3n4);P2q2;— ptqt— ;Prqr 
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Similarly, by fp^Ojl and (j02]l . 

jja6,ZiZ2,( n i n 2"3|;Pi,|n4)qi2gi3;P2q2i— ;Prqr 

_ rpo6;Zii2,(nin2n3|;pi,|ri4)?i2gi3;p2q2;— ;Prqr + ^W2;at>,(niri2n3|;pi,|«.4)<2i2<ji3;P2q2;---;Prqr 
_ rpgi2gi3;'i'2,(™ir»2ri3l;Pi,|n4)afe;P2q2;---;Prqr _j_ rpgi2gi3;a6,(nin2n3];pi,|n4)ii22;P2q2;-";Prqr 

By the above, equation (14.231) becomes 

(m + 2)T <?12 ' ?13 ' a ''' PllPl2( - ni ''' 1 ' 2 '' n2r '' 3n4 ^ ;P2q2; "' ;p, ' qr 

r 

+ V* r p9i29i3iab,PtiPt2(ni|;2ii2,|n27i3n4);piqt;p2q2i---ptqf"iPrqr 

h (4.24) 
_l_ ^912913 ;W2,(nin2n3|;pi,|n4)af>;p2q2;---;prqr 

+ rp5i2gi3;a6,("in.2n3|;pi,|n4)h^2;P2q2;---;Prqr _ g 

We have thus reduced (I4.23P into an equation for the derivatives of a single component 
T 912913 of the source form T. 

Now assume that equations (I4.24p admitted a non-trivial solution in the derivatives 
of T 912913 . Choose a term T 9l29l3;af '' piipi2(ni|; ' li2 ' |ri2n3?14);p2q2; - ;p '- q '' with maximal absolute 
value amongst the symmetrized derivatives. Keeping in mind that r < m - 2 by the 
induction assumption, equation (14.241) implies that 

| T ft 3 ff U ;a6 >Pll p la (ni|;l 1 I 3 ,|n 9 nsru);p a q 3 i™;Prqr| 
1 



gi2gi3;a6,PtiPt2(ni|;'i'2,|n2"3"4);piqt;P2q2;--ptqf;Prqr 

t=2 
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m + 2 1 

+ | T ffi ag is;iii a ,(n 1 n a n3| ;P1 ,|n4)a6;p a w..;Prqr| + | r j i ai2gi3;afe,(nin2n3|;pi,|n.4)^i^2;p2q2;--;prqr | j 

< m ~ 1 |rpg 2 g 3 ;a6.p}p 2 (ni[;M2,|ra2n3n 4 );p2q2;---;Prqr| 

~ m + 2 1 '' 
Thus necessarily 

jjafe,cd,e/;ii/ 2 ,nin2n3™4;Piqi;---;Pr-iqr-i _ rjie/;afc,cd(ni|;M2 J«2™3™4);Piqi;---;Pr-iqr-i _ g ^ 

for all a, b, c, d, e, f, h, l2,ni, . . . , n^, pi, qi, . . . , p r _i, q r -i- Hence (I4.2ip reduces to 

^(ni|jjafe,/i/2,|n,2n3n-4);Piqi;---;Prqr _ g 

which immediately implies that H° 6 > /l/ 2> ra i rt 2™.3,piqir--,Prqr _ g^ com pi e ^i n g the induction 
step. 

In conclusion, the third order components }j^> cd >i kl Q f the Helmholtz operator asso- 
ciated with T vanish identically. By Corollary [TT] the Helmholtz operator Hx must also 
vanish identically, and hence T is locally variational. \qed\ 
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5 Discussion 



In this paper we prove that a system of third order natural, divergence- free differen- 
tial equations for the components of the metric tensor can always be written as the 
Euler-Lagrange expression of some Lagrangian function. Moreover, by the solution of 
the equivariant inverse problem of the calculus of variations for metric field theories 
presented in [5], the Lagrangian can be chosen to be natural when the dimension of the 
underlying space is m = 0, 1, 2 mod 4, and, in dimensions m = 3 mod 4, the system can 
be written as a sum of the Euler-Lagrange expression of a natural Lagrangian and a 
generalized Cotton tensor. 

Extending our result to fourth order operators either by showing that natural, diver- 
gence free systems are necessarily locally variational or by finding non- variational ex- 
amples of such systems remains a challenging open problem. In contrast, for vector field 
theories |8| and, in more generality, for Yang-Mills theories |21j, non- variational third 
order source forms with the prescribed symmetries and conservation laws can be de- 
rived by the way of natural constructions utilizing the intrinsic geometry of the problem. 
Note that unlike in the situation of the present paper, proving that a 4th order source 
form T fulfills the highest order non- vacuous Helmholtz conditions H^' ccMli2W4 = will 
not be sufficient to guarantee that T is locally variational. 

A significant generalization of the present work would be to enlarge the symmetry 
group to include conformal transformations of the metric with the resultant condition 
that the source form be trace free. The problem closely bears on conformal gravity, 
and, as exemplified by the Bach equations pU], will require the analysis of fourth order 
systems. A solution for the equivariant inverse problem of the calculus of variations in 
this situation would also be of substantial independent interest. Additional on physi- 
cal grounds intriguing problems would be Takens' question for combined field theories 
involving metric, such as Einstein- Yang-Mills equations in both the abelian and non- 
abelian cases. 
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